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ABSTRACT: Pork from Jeju black pig (population J) and Berkshire (population B) has a unique market share in Korea because of 
their high meat quality. Due to the high demand of this pork, traceability of the pork to its origin is becoming an important part of the 
consumer demand. To examine the feasibility of such a system, we aim to provide basic genetic information of the two black pig 
populations and assess the possibility of genetically distinguishing between the two breeds. Muscle samples were collected from 
slaughter houses in Jeju Island and Namwon, Chonbuk province, Korea, for populations J and B, respectively. In total 800 Jeju black 
pigs and 351 Berkshires were genotyped at thirteen microsatellite (MS) markers. Analyses on the genetic diversity of the two 
populations were carried out in the programs MS toolkit and FSTAT. The population structure of the two breeds was determined by a 
Bayesian clustering method implemented in structure and by a phylogenetic analysis in Phylip. Population J exhibited higher mean 
number of alleles, expected heterozygosity and observed heterozygosity value, and polymorphism information content, compared to 
population B. The F ls values of population J and population B were 0.03 and -0.005, respectively, indicating that little or no inbreeding 
has occurred. In addition, genetic structure analysis revealed the possibility of gene flow from population B to population J. The 
expected probability of identify value of the 13 MS markers was 9.87x10 14 in population J, 3.17xl0~ 9 in population B, and 1.03x10 12 
in the two populations. The results of this study are useful in distinguishing between the two black pig breeds and can be used as a 
foundation for further development of DNA markers. (Key Words: Black Pig Populations, Genetic Traceability, Heterozygosity, 
Microsatellite Markers, Probability of Identity) 



INTRODUCTION 

In Korea, fresh pork is popular for use in a wide variety 
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of dishes. With this arises a demand for high quality fresh 
pork with high marbling. In Japan, the retail price of 
purebred Berkshire meat is 50% higher than that of regular 
commercial pork meat (three-way crossbred pigs) (Suzki et 
al., 2003). Similarly in Korea, black pork, mostly from 
Berkshire and Jeju Black pig, is sold at a higher price than 
regular commercial pork. 

Previous studies of black pig meat revealed a higher 
glutamic acid content and higher red muscle fiber content 
than other commercial breeds. The level of glutamic acid is 
known to be a major factor in flavor (Moon et al., 2004) 
while the red muscle fiber content is thought to confer 
tenderness to the meat (Song et al., 2010). In addition, it 
was reported that black pig meat, especially from Berkshire, 
possesses high marbling, drip loss and water holding 
capacity (Lee et al., 2011). Due to these characteristics, the 
black pig meat is in high demand in Korea. 

Recently, the Korean public has started to pay attention 
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to the origin of the food it consumes. To provide the 
consumers with the necessary information, the potential of 
using DNA-based methods for the identification of meat 
has been extensively explored (Dalvit et al., 2007; 2008). 
Several types of DNA markers have been tested for use in 
traceability systems in livestocks, including mitochondrial 
DNAs (Alves et al., 2009), microsatellites (Kim et al., 
2007), and SNPs (Goffaux et al., 2005; Negrini et al., 2008; 
Shim et al., 2010). A traceability system based on DNA 
technology has already been enforced for Korean beef since 
2009 in Korea. To consider expanding this use of genetic 
identification system to pork, the Korean government has 
formed a committee to consider suitable methods for 
implementing a traceability system. In addition, multiple 
providers supply the market with black pig meat in Korea 
and the use of a genetic system to distinguish between 
different stocks will also enhance the brand power of the 
suppliers. 

With the goal of developing a genetic tool for the 
identification of pig meat, the objective of this study was to 
identify markers with sufficient genetic information to 
discriminate between the two main black pig populations in 
Korea. A total of 13 microsatellite (MS) markers were 
tested in Jeju black pig and Berkshire, and levels of genetic 
diversity in each breed were estimated for future use in 
breeding. 

MATERIALS AND METHODS 

Sample collection and DNA extraction 

Muscle samples were collected between September of 
2012 and April of 2013 from two slaughter houses; one in 
Jeju Island for population J (n = 800) and the in Namwon, 
Korea, for population B (n = 351). Population J is known to 
have a mixed genetic background with possible cross- 
breeding with either Duroc or Berkshire and Korean native 
pig. Population B consists of pure Berkshire breed, which 
was established in Korea since 2004. 

Total genomic DNA was extracted using the salting-out 
procedures (Miller et al., 1988). DNA quality and quantity 
were measured with UV-spectrophotometer. 

Microsatellite marker genotyping 

In this study, we selected 13 swine microsatellite 
markers to analyze the genetic polymorphisms of the two 
pig populations (Table 1). Among 13 markers, 10 markers 
were recommended by ISAG/FAO (2011) and 3 markers, 
i.e., S0225, SW951 and SW787, which were recommended 
by Roslin Institute (University of Edinburg, UK) for pig 
genetic diversity studies. Also, these markers were 
previously used for individual identification and parentage 
analysis in commercial pigs in Korea (Lim et al., 2009). 
Polymerase chain reaction (PCR) was performed in 15 uL 



Table 1. List of microsatellite markers used in this study 



Locus 


Chromosome 


Dye 


Size range (bp) 


S0155 


1 


NED 


145-168 


SW72 


3 


NED 


95-134 


S0005 


5 


VIC 


200-257 


SW122 


6 


VIC 


116-138 


SW632 


7 


PET 


151-190 


S0225 


8 


NED 


-t f~tf\ o r\f\ 

172-200 


SW951 


10 


FAM 


123-145 


S0090 


12 


FAM 


200-260 


SW857 


14 


VIC 


142-168 


SW936 


15 


FAM 


80-120 


S0026 


16 


VIC 


90-113 


SW24 


17 


PET 


91-129 


SW787 


18 


FAM 


150-180 



volume containing of 1 uL of 100 ng of genomic DNA, 0.5 
|iL of 10 pmole of each primer, 1.5 |iL of lOxPCR buffer 
(10 mM Tris-HCl, 50 mM KC1, 1.5 mM MgCl 2 , pH 8.3), 
1.5 uE dNTP (2.5 mM), and 0.2 uE Taq polymerase (10 
unit/ uL) (GENETBio, Daejeon, Korea). Thermal cycling 
conditions in GeneAmp PCR System 9700 (ABI thermal 
cycler, Foster city, CA, USA.) included an initial 
denaturation for 5 min at 94°C, followed by 35 cycles of 30 
s at 94°C, 30 s at the annealing temperature, 1 min at 72°C, 
and a final extension step of 10 min at 72°C. The amplified 
products were visualized on 2% agarose gels stained with 
ethidium bromide (0.5 ug/mL). Genotyping was performed 
on ABI 3130 Genetic Analyzer (Applied Biosystems, Foster 
city, CA, USA). 

Statistical analyses 

The number of alleles, observed heterozygosity (H 0 b), 
expected total heterozygosity (H Ex ), and polymorphic 
information content (PIC) were calculated using MS toolkit 
software (Park, 2000). Pairwise F ST and inbreeding 
coefficients (F IS and F IT ) were calculated using the program 
FSTAT (Goudet, 2001). The clustering method (Saitou and 
Nei, 1987) was used to construct a phylogenetic tree based 
on the genetic distance matrix in Phylip package version 
3.67 (Felsenstein, 2007). To assign individuals to 
populations using genotype data, a Bayesian approach was 
used in the program structure (Pritchard et al., 2000). The 
expected probability of identity (PI) value among genotypes 
of random individuals was calculated using API-CALC 
version ver 1.0 (Ayres and Overall, 2004). 

RESULTS AND DISCUSSION 

Genetic diversity of the two black pig populations 

We genotyped 1,151 black pig samples, 800 Jeju black 
pigs (J) and 351 Berkshires (B), with 13 MS markers. These 
markers were recommended by ISAG for both breed and 
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Table 2. Number of samples, mean number of alleles (MNA) per locus, expected heterozygosity (H Ex ) at Hardy- Weinberg equilibrium, 
observed heterozygosity (H 0 b), polymorphism information content (PIC) and inbreeding coefficient (F IS ) of the two black pig 
populations 



Population 


No. of samples 


MNA 


Hex 


H 0 b 


PIC 


F,s 


J 


800 


8.23 


0.706 


0.681 


0.668 


0.030 


B 


351 


5.08 


0.541 


0.544 


0.486 


-0.005 



J, Jeju black pigs; B, Berkshire. 



individual identification in pigs. Furthermore, the markers 
have been validated by a previous study on breed 
identification and origin tracking for pork meat in Korea 
(Lim et al., 2009). 

The expected H Ex , H 0 b, and PIC are important 
parameters for the selection of markers for pig breed 
identification. In general, mean number of alleles per locus 
(MNA), heterozygosities (H Ex and H 0 b), and PIC value of 
population J were higher than those of population B (Table 
2). The MNA of population J and B was 8.23 and 5.08, 
respectively. H Ex of population J was 0.706, whereas that of 
population B was 0.541. Population J showed a higher PIC 
value than population B (0.668 vs 0.486). 

Genetic makers for identification of breed need to have 
sufficient genetic diversity and power. A previous study has 
indicated that for successful individual identification, PIC 
and H Ex values need to be higher than 0.5 and 0.6, 
respectively (Botstein et al., 1980). 

H Ex of the two pig populations combined ranged from 
0.416 (S0225) to 0.855 (S005), and PIC value ranged from 
0.338 (S0225) to 0.840 (S005) for individual markers. 
Mean H Ex and PIC were 0.688 and 0.650, respectively 
(Table 3). A few markers, S0225, SW951, SW026, appeared 
to be less polymorphic than others. On the contrary, some 
markers, such as S005, SW857, SW122, and SW632, 
showed higher genetic diversity and information making 



them especially useful for traceability tests. In population J, 
H Ex values ranged from 0.484 (SW951) to 0.884 (S005), 
while in population B, they ranged from 0.146 (SW0225) to 
0.772 (SW24) (Table 4). The highest PIC value was seen in 
marker S005 (0.873) for population J, while for population 
B, it was marker SW24 (0.743) (Table 4). Marker S0225 
showed the lowest PIC value for both population J and P 
with 0.447 and 0.136, respectively. These analyses 
demonstrated that 10 out of 13 MS markers tested may 
provide a reasonable amount of genetic diversity and 
information to be suitable for use in breed identification. 

Analyses of Wright's F-statistics were carried out using 
FASTA (version 2.9.3) and the results are presented in Table 
3. The global heterozygosity deficit between two black pig 
populations was estimated by F IT using the 13 MS markers. 
The mean F IT and F ST among the markers were 0.126 and 
0.194, respectively. F lT value was the highest for SW122 
(0.191), while it was the lowest for SW0026 (0.022). F ST 
value of SW122 was the highest (0.175) and that of 
SW0026 was the lowest (0.026) (Table 3). 

The inbreeding coefficient, F IS , reflects the degree of 
inbreeding in a population by measuring the loss of 
heterozygosity within a population. F is ranged from -0.054 
(SW857) to 0.120 (S0005) and the mean F IS of the 13 
markers was 0.025 (Table 3). Table 4 shows the F\ S values 
for population J (0.03) and population B (0.005), which 



Table 3. Measures of genetic diversity for the 13 microsatellite markers 


Locus 


Hex 


Hob 


PIC 


Fyt 


^ST 




S0155 


0.681 


0.690 


0.649 


0.027 


0.068 


-0.044 


SW72 


0.662 


0.597 


0.609 


0.129 


0.059 


0.075 


S0005 


0.855 


0.736 


0.840 


0.165 


0.051 


0.120 


SW122 


0.792 


0.713 


0.760 


0.191 


0.175 


0.019 


SW632 


0.793 


0.763 


0.760 


0.090 


0.093 


-0.002 


S0225 


0.416 


0.402 


0.382 


0.105 


0.129 


-0.027 


SW951 


0.474 


0.425 


0.448 


0.139 


0.068 


0.076 


S0090 


0.731 


0.676 


0.684 


0.167 


0.161 


0.006 


SW857 


0.748 


0.716 


0.701 


0.151 


0.194 


-0.054 


SW936 


0.706 


0.670 


0.655 


0.087 


0.065 


0.024 


S0026 


0.578 


0.573 


0.520 


0.022 


0.026 


-0.003 


SW24 


0.816 


0.710 


0.792 


0.153 


0.046 


0.112 


SW787 


0.694 


0.650 


0.652 


0.155 


0.169 


-0.017 


Mean 


0.688 


0.640 


0.650 


0.126 


0.103 


0.025 



Hex, expected heterozygosity; H 0 b, observed heterozygosity; PIC, polymorphism information content; Frr, measures of the deviation from Hardy- 
Weinberg in the total population; F ST , measures of the genetic differentiation over subpopulations; Fi S , measures of the deviation from Hardy- Weinberg 
proportions within subpopulations. 
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Table 4. Statistical analysis of heterozygosity (H Ex and H 
populations using 13 microsatellite markers 


ow> riL ana 


gene differentiation within population (F IS ) of two black pig 


Locus 








n Ob 




PIC 


^IS 




J 


B 


J 


D 
D 


J 


B 


J 


B 


S0155 


0.692 


0.585 


0.700 


U.OOJ 


0.661 


0.512 


-0.011 


-0.136 


SW72 


0.634 


0.669 


0.550 


U. /Vj 


0.591 


0.595 


0.133 


-0.051 


S0005 


0.884 


0.718 


0.759 


U.Ool 


0.873 


0.666 


0.142 


0.051 


SW122 


0.799 


0.559 


0.793 


\j.jZy 


0.769 


0.526 


0.009 


0.054 


SW632 


0.783 


0.711 


0.769 


n 7/1 Q 
U. /4y 


0.747 


0.657 


0.018 


-0.055 


S0225 


0.499 


0.146 


0.514 


U. 140 


0.447 


0.136 


-0.030 


-0.002 


SW951 


0.484 


0.404 


0.469 


0.325 


0.447 


0.356 


0.031 


0.197 


S0090 


0.727 


0.552 


0.724 


0.544 


0.679 


0.467 


0.004 


0.015 


SW857 


0.736 


0.550 


0.753 


0.631 


0.686 


0.496 


-0.023 


-0.148 


SW936 


0.742 


0.556 


0.740 


0.509 


0.699 


0.457 


0.003 


0.086 


S0026 


0.619 


0.462 


0.608 


0.494 


0.562 


0.404 


0.018 


-0.069 


SW24 


0.812 


0.772 


0.700 


0.734 


0.785 


0.743 


0.138 


0.049 


SW787 


0.767 


0.347 


0.779 


0.357 


0.731 


0.302 


-0.015 


-0.029 



J, Jeju black pigs; B, Berkshire. 



suggests that little or no inbreeding has occurred within the 
two populations. This is most likely due to the fact that both 
populations have been under managed breeding systems. 

Phylogenetic and structure analyses 

An individual-based neighbor-joining dendrogram was 
built from the estimates of allele-sharing distances among 



the 400 individuals which consisted of 200 Jeju black pigs 
and 200 Berkshires (Figure 1). Based on the equation of Nei 
et al. (1983), a phylogenetic tree was estimated by the 
distribution of allele sharing by genetic distance (F ST ). The 
two black pig populations were clearly distinguishable by 
MS markers. Population J showed a large genetic distance 
from population B. This indicates that although the two 



J535 




Figure 1. A neighbor joining dendrogram was constructed from allele-sharing distances among 400 individuals in the two pig 
populations. J: Jeju black pigs, B: Berkshire pigs. 



930 



Oh et al. (2014) Asian Australas. J. Anim. Sci. 27:926-931 



populations share a coat color, they are genetically distinct 
from each other. We observed that the distance between the 
branches in population J was greater than that of population 
B. This suggests that population J should have a greater 
variation than population B. 

We estimated the expected PI values of the 13 MS 
markers in random individuals in the two black pig 
populations (Table 5). The expected PI from was 9.87xl(T 14 
and 3.17xl(T 9 for population J and population B, 
respectively. Kim et al. (2007) reported that analysis among 
4 pig breeds (Duroc, Landrace, Yorkshire and Korean native 
pig) with 10 MS markers resulted in an expected probability 
of identify value of 0.36xl0~ 9 . In our study, the expected 
probability of identify value was estimated to be 3.87x10"" 
with 10 markers. In another study of 400 pigs with a three- 
way cross (LandracexLarge WhitexDuroc), the expected 
probability of identify value was estimated to be 2.47 xl0~ 18 
(Lim et al., 2009). Therefore, the results of this study are 
reliable when compared to previous studies. In summary, 
the expected PI value was 1.03xl0~ 12 with 13 MS markers 
in 1,151 black pigs (Table 5). 

The result of the structure analysis when K = 2 is shown 
in Figure 2. Each breed grouped into its own cluster; 
however, there were signs of gene flow from population B 
into population J. The possibility of cross breeding between 
the two populations is supported by the breeding history 
(Moon, 2004). 



Table 5. The expected PI values among genotypes of random 
individuals for discriminating between the two pig populations 



INO. Ul IlldlKcIS 


T 
J 


R 
O 


lUlal 


AA 1 


Z. 1 1X1U 


i ^/w 1 n - 2 

/.JOXIU 


J.ZjXIU 


1V1Z 


1 i ^ 1 n - ^ 
1 . 1 jX iu 


/ AjXLV 


1 .o4X!U 




/.UjXIU 


q i Aw i 
o. 1UX iu 


1 .j 1X1U 


AAA 


z 1 1 w i n~6 
J. 1 1X1U 


i m v 1 
1 .uzxiu 


v.^yxiu 


M5 


4.07x1 0" 7 


1.94xl0" 5 


9.08xl0" 7 


M6 


3.96xl0" 8 


4.03xl(T 6 


9.96xl0" 8 


M7 


4.11xl0" 9 


8.49xl(T 7 


1.35xl0" 8 


M8 


4.38xl0" 10 


2.42xl0" 7 


1.83xl0" 9 


M9 


4.87xl0 _il 


6.10xl0 -8 


2.32xl0" 10 


M10 


7.56xl0" 12 


1.99xl0" 8 


3.87xl0" 11 


Mil 


1.39xl0" 12 


9.63xl0" 9 


8.95xl0" 12 


M12 


3.71xl0" 13 


4.43xl0" 9 


2.77xl0" 12 


M13 


9.87xl0" 14 


3.17xl0" 9 


1.03xl0" 12 



PI, probability of identity; J, Jeju black pigs; B, Berkshire. 



conventional traceability system based on meat cut labeling. 
A genetic traceability system will be useful for identifying 
black pig meat products in the market to prevent fraud and 
mislabeling. However, trade-off between economic viability 
and sufficient identification power need to be considered for 
the selection of a subset of MS markers for a practical 
traceability system. 
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